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Abstract 



ON ■ 
O ; 

Although there is no known meaningful notion of the energy density of the grav- 
\ itational field in general relativity, a few notions of quasi-local energy of gravity 

associated to extended but finite domains have been proposed. In this paper, the 
notions of quasi-local energy are studied in the framework of loop quantum gravity, 
in order to see whether these notions can be carried out at quantum level. Two basic 
quasi-local geometric quantities are quantized, which lead to well-defined operators 
' ^j' in the kinematical Hilbert space of loop quantum gravity. We then use them as 

Q" 1 basic building blocks to construct different versions of quasi-local energy operators. 

5_i ' The operators corresponding to Brown- York energy, Liu-Yau energy, Hawking en- 

, ergy, and Geroch energy are obtained respectively. The virtue of the Geroch energy 

operator is beneficial for us to derive a rather general entropy-area relation and thus 
| a holographic principle from loop quantum gravity. 

lO : PACS number(s): 04.60.Pp, 04.20. Cv 

in ■ 

<N : 1 Introduction 

oo : 

O ■ It is well known that there are inherent difficulties in defining energy in general relativity 
J> ! (GR), essentially owing to its non-localizability. By now there is no known meaningful 

notion of the energy density of gravitational field in GR. Globally, for spacetimes which 
are asymptotically flat, there are well-defined notions for the total energy, given by the 
Bondi and ADM expressions integrated over spheres at null infinity and spatial infinity. 
These global notions are directly related to quantities that can be measured physically by 
distant observers. However, finding an appropriate notion of energy-momentum would be 
important from the point of view of applications as will. For example, the correct, ultimate 
formulation of black hole thermodynamics should probably be based on quasi-local defined 
internal energy, entropy, angular momentum etc. So far, considerable efforts have been put 
in to formulate a satisfactory definition of quasi-local energy (QLE) (see [1] for a review). 
In this paper, a few expressions of quasi-local energy are quantized in the framework of loop 
quantum gravity (LQG) (see [2-5] for reviews). Our purpose is in two folds. Firstly, we 
want to check whether the quasi-local notions of gravitational energy can be carried out at 
quantum level. Secondly, we wish to use these notions of quantum gravitational energy to 
study the relation between quantum gravity and gravitational thermodynamics. A similar 



*e-mail address: yangksong@gmail.com 
' e-mail address: mayg@bnu.edu.cn 



1 



effect was made in [6] to quantize the Hamiltonian surface term of a bounded spatial region. 
The candidates of QLE that we are considering include the Brown- York energy [7], Liu-Yau 
energy [8], Hawking energy [9], and Geroch energy [10]. All of these energy expressions 
are constituted by two basic quasi-local quantities representing the extrinsic curvatures of 
an two-sphere in a spatial slice and a timelike slice respectively. Thus the key task is to 
quantize these two basic building blocks. 

In Section 2, the kinematical framework of LQG is briefly introduced. We then con- 
struct in section 3 two basic operators in the kinematical Hilbert space of LQG, which 
corresponding to the two basic building blocks of quasi-local energies. In section 4, the two 
basic operators are employed to construct different versions of QLE operators. In section 
5, the Geroch energy operator is used to derive a rather general entropy-area relation and 
thus a holographic principle from LQG. 



2 Elements of LQG 

The Hamiltonian formalism of GR is formulated on a 4-dimensional manifold M = R x 
E, where E represents a 3-dimensional manifold with arbitrary topology. Introducing 
Ashtekar-Barbero variables [11,12], GR can be casted into an SU(2) connection dynamical 
theory. The phase space consists of canonical pairs (A^, Ef) of fields on E, where A l a is a 
connection 1-form which takes values in the Lie algebra su(2), and Ef is a vector density 
of weight 1 which takes value in the dual of su(2). Here a, b, c... are abstract spatial indices 
and i, j, k... = 1, 2, 3 are internal su(2)-indices. The density-weighted triad Ef is related to 
the co-triad e l a by the relation Ef = \ e ahc eiji.e b e^ sgn(det(ejj), where e abc is the naturally 
defined levi-civita density and sgn(det(ejj) denotes the sign of det(e^). The 3-metric on E 
is related to the co-triad by q a b = e l a e J b 5ij. The only non-trivial Poisson bracket is given by 

{Al(x),E b (y)} = ^5 b a 5}5 3 (x,y), (2.1) 

where k = 8ttG (G denotes Newton's constant) and (3 is the Barbero-Immirzi parameter. 
There are three first-class constraints in this Hamiltonian formalism of gravity: 



Gi = V a Ef = d a E? + e ijh A{E\ 



hi 

V a = Fl h El 

rpaj rpbk 

H = —^=== \e iik F! b - (1 + f3 2 )2Kl Kt)) , (2.2) 
2«y|det(g)| L 3 ab 1 ' [a bVi 1 ' 

where T> a denotes the covariant derivative defined by the connection A l a , F* b := d a A\ — 
dbA l a + e l jkA J a A b is the curvature of A l a , and K l a is the extrinsic curvature of E. 

One element of LQG is the notion of graphs embedded in E. By 7 we denote a closed, 
piecewise analytic graph. The set of edges of 7 is denoted by ^(7) and the set of vertices 
of 7 by V(7). For an oriented edge e of 7, its beginning point is denoted by b(e) and its 
final point by /(e). To construct quantum kinematics, one has to extend the configuration 
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space A of smooth connections to the space A of distributional connections. Through 
projective techniques, A is equipped with a natural, faithful, 'induced' measure /x G , called 
Ashtekar-Isham-Lewandowski measure [13,14]. In certain sense, this measure is the unique 
diffeomorphism-invariant measure on A [15]. The kinematical Hilbert space then reads 
7ikin = L 2 (A,d/i ). The so-called spin- network basis T 7i j, m ,n provide an orthonormal basis 
for H km [2]. 

One of successes of LQG is the rigorous construction of spatial geometrical operators, 
such as the area, the volume and the length operators in 7Ykm [16-18]. Moreover it turns 
out that these geometrical operators have a discrete spectrum. The same conclusion are 
also tenable in the internal gauge invariant Hilbert space 7i = L 2 (A/G, d/j, ). 

3 Two basic operators for QLE 

Most of the quasi-local energy expressions appeared so far involve two quasi-local quantities 
defined by the integrals of two extrinsic scalar curvatures of some spatial 2-surface [1]. In 
this section, we will construct two well-defined basic operators corresponding to the two 
quasi-local quantities in the kinematic Hilbert space 7Ykin of LQG, which will be used as 
basic building blocks to construct different versions of quasi-local energy operators in the 
next section. To this aim, we have to first re-express the two quasi-local quantities in terms 
of real connection variable or its conjugate. Then we regulate the classical expressions in 
order to get quantities with quantum analogues. It turns out that in the regularization 
procedure, as the regularization of the Hamiltonian constraint, we need to triangulate the 
3-d spatial manifold £ in adaption to a graph, which comes from the cylindrical function 
in Hkin that is going to be acted by the constructed operators. 

3.1 QLE-like operator 

Let S be a 2-d surface with two-sphere topology in the 3-d spatial manifold £ and a ab be 
the induced metric on S of metric q ab on S. For simplicity, we choose adapted coordinates 
{x 1 , x 2 , x 3 } in £ such that S is given by x 3 = 0, and x x ,x 2 parameterize S. The QLE-like 
observable is defined as 

E Q , k (S) := -- [ d 2 x^Jdet(a) k, (3.1) 
K Js 

where k is the extrinsic scalar curvature of k a b of S corresponding to the unit normal n a 
in £. In order to quantize the expression, we need first to express it in terms of the real 
Ashtekar variables. The extrinsic curvature tensor k ab of S corresponding to n a reads 

k ab = a a c D c n b , (3.2) 

where D a is the derivative operator on S compatible with q ab , i.e., D a q bc = 0. So the 
extrinsic scalar curvature k of S is 

k =D a n a . (3.3) 
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From (13. 3p . we have the following identity in the adapted coordinates. 

/ E a E b \ 1 f E a E b \ 1 ( E 3 E a 

k = D A^r\ nb ) =^t=tt^ \-7T=h n » =-^=tt 9 - NI=M- ( 3 - 4 ) 

\det(q) J y/det(q) \y/det(q) J y/det(q) \y/det(cr) 

To regularize Eq^(S), let e be a small number and (x, y) be the (smoothed out) character- 
istic function such that \im € ^ x 3 (x,y)/e 3 = 5 3 (x,y). The volume of the cube as measured 
by q ab is given by V(x, e) := / d 3 yx 3 (x, y)y/det(q)(y) such that lim e ^ = y/det(q)(x). 
We choose again adapted coordinates {y 1 , y 2 , y 3 } in £ such that S is given by y 3 = and 
each 2-d surface S y a of the family given by y 3 = constant is parameterized by y x ,y 2 . We 
denote also the induced metric of q a b in S y 3 by <r a b. We then have the following identity 
by Eqs. (ED) and (E3D 

QM ' kJ s ^/detM 

= - lim - / rf 2 x v/delRx) / */ 3 / rfy 1 ^ 2 <9 a (j%SL-(y) 

t-^KJs J Js e 3 \ v/det(a) 



x / 



e 
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Vdet( g )(u) 



lim- / d 2 x^t(a)(x) [ dy 3 [ dy l dy 2 x 3 e (x,y) d a ( ^2= (y) 
K Js J Js., 3 \ \/det a 



v/det( 



fi...*t„ ..n [det(e|)] («) /■ * , [det(e^)] (w) 



x / d s u X i(x,u) ) b,i \ ' / d 3 w X i(x,w) 
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Hm- / rf 2 xv^et(^)(x) / dy 3 [ dy l dy 2 [d aX 3 (x,y)] [ -5^L(y) 
A., ,3/„ .a t det (4)] («) f * ...x [det(e|)] H 



where we have inserted the identity 1 = [det(e^)] 2 /[^/det(g) ] 2 in the third step, and 
performed an integration by parts in the last step . Let x\> {Vi z ) is the 2-d characteris- 
tic function of a coordinate box with center y and coordinate area e' 2 and Ar(z, e') := 
fc d 2 wx\i (z,w)y/ det (a) (it;) , satisfying rim e /_> Ar , 2 ' e ^ = */ det(cr)(,z), is the area of the 
box as measured by a a b. Then we have 

Ef i 1 2 3\ r /" „r2 Xe'(Z/ ' 2/ 2 ; 2 i ^ I 1 2 3\ 



v/d^) ^0^3 e ' 2 v /5 ^) 



lim /• d2z x 2 Ay\y 2 ;z\z 2 )E?(z\z 2 ,y 3 ) 

™ Js y3 Z Ar(z\z 2 ,y 3 ;e>) 
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Thus we can rewrite Eq. f l3.5j) as 

En.k{S) = lim - / d 2 xy/det {a)(x) 

e,e'^0 K Js 



x / c? 3 2/ 



> Xe'iv 1 ^ 2 ; z 1 , z^Efiz 1 , z 2 ,y 3 



d z z 



S y 3 



Ar(z 1 , z 2 ,y 3 ; e') 



[det(ej)] (u) 



Recall the following classical identities 



d 3 wx 3 e (x,w) 



[d a &x,y)] Efty) 
[det(4)] (w) 



(3.6) 



rf 3 w [det(e l 6 )] (u) = - / e ijfc e 4 A e J ' A e fc = — / Tr(e A e A 



3! 



and 



eU M ) = -K(n),y(n,6)}, 



(3.7) 



(3.8) 



where e^. = — ~ tr(r i r ? r fc ), e = eVj/2, here 7$ = — icr^ (<t» is the Pauli matrix) is the 
generator of sw(2) obeying [r^Tj] = 1e^ k T k . We can rewrite Eq. f!3.6p as 



-4 /2\3 



3! Vk 
x / d 3 y 



, X ?,(^VV 2 )^(*V 2 ,</ 3 ) 



S H 3 



Ar(z 1 , z 2 , y 3 ; e') 



[9 aX ?(x,y)]£?(y) 



v v 7 ^) \/%) 

y E AeV ^ aA>^) Vvm 



2 1A r 

— - lim / d 2 x\/det(a)(x) 



x / d 3 y 



[Sax?(*,y)]^?(y) 



Ar(z 1 , z 2 , y 3 ; e') 

x / x 3 (x, u)Tr v 7 ^, e) } A {A(u), y/V(u, 7) } A {A(u), y/V(u, e) } 

x / x 3 (x, w)Tr [{A(w), y/V(w, e)} A {A(w), y/V(w, e) } A {A(w), y/V(w, 7) } 
--: \un o E^ k (S), (3.9) 
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where we have used {■, ^/V(u, e)} = {•, V(u, e)} / {2^/V{u, e) ). Since a typical state 
/ 7 G TYkin is some cylindrical function over a graph 7 in £, as in the construction of 
the Hamiltonian constraint operator [2,19], we triangulate £ in adaption to 7 as follows. 
At every vertex v G V(7) we choose a triple (ej, ej, e#) of edges of 7 and a tetrahedron 
A%,ej,ej,e K based at w which is spanned by segments si, sj, sk of the triple. Each segment 
Si is given by the part with the curve parameter t 1 G [0, e] of the corresponding edge e^t 1 ). 
The holonomy of the connection along a segment Si reads 

h Sl (A) =I 2 + es a j(0)Al(v) n /2 + 0(e 2 ), (3.10) 

and for one segment sj, we have 

{A(u), VVM} « es?(0){A,(«), VVM) (3.11) 

up to 0(e 2 ). Hence for each A* , we have 

'A(u), y/V(u,e) } A [A(u), y/V(u,e) } A y/Vfae) }) 

^e( S/Sj ^)6 7JX Tr(/ isj(A) {/i- 1 ^, ^(«(A),e) } h Sj{A) [h^ (A) , ^/V(v(A), e) } 



Tr 



X /W(A) {^(A)> V^K%) } ): 



where e(s/SjSft:) := sgn(det(sfSjs^)(i;)) takes the values +1,-1,0 if the tangents of the 
three segments si,sj, sk at v (in that sequence) form a matrix of positive, negative or 
vanishing determinant. Then the integration over £ can be split as follows [2,19]: 



+ £ L = /_+ £ 



S ^ ..nsf^JUZ,, JU. 



1 veV(j) u -i,v u i v£V(~j) v ' b(e I )nb(e J )nb(e K )=v 



+ 



7 VGV("/) y ' b(e I )nb(ej)nb(e K )=V L '" i 7,t>,ej,ej,e Jf U; f,v,ej ,ej,e K _ 



(3.12) 



Here we have first decomposed E into a region £7* not containing the vertices of 7 and the 
regions around the vertices. Then choose a triple (ej, ej, ck) of edges outgoing from 
v and decompose Ut n . into the region U%,..^.^.„ covered by the tetrahedron Al„ _ , „ T<> „ 
spanned by ej, ej, and its 8 mirror images and the rest ei ej ex not containing v. 
Note that the integral over U* ei ej eK classically converges to 8 times the integral over the 
original single tetrahedron Af, v e/ ej eK as we shrink the tetrahedron to zero. We average 
over all such triples (ej, ej, ex) and divide by the number of possible choices of triples for a 

vertex v with n(v) edges, E{v) = ("^V We can now decompose the u and w integration 

over £ in Eq. (13.91) according to Eq. (I3.12p . To quantize the two integration, we replace 
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Poisson brackets by commutators times l/(ih), holonomies by multiplication and V by the 
volume operator V, which acts on a function cylindrical over a graph 7 as follows [17]: 



V(R)f y = (W?) 3/2 E 



\ 



3! -8 



E e(e,e',e'% jk XlXi,X k e „ 



ene'r\e"=v 



(3.13) 



Because the non-vanishing contributions of V acting on a cylindrical function / 7 come 
from the vertices v G ^(7), only the integration over the tetrahedra A„ vej ejeK needs to 
be considered. Hence the resulted operator corresponding the last two integrals in Eq. 
(I3.9P acts on a cylindrical function as 



36z 6 fr 6 



E * 



»'ev(7) 



- ^ e(s/SjSx)e /J ^Tr(/i S/( 



u(A)=t/ 



s/(A)' 



V>',e) 



x ft, 



*j(A) 



sk(A) 



X 



E xfc^) 



«"6V(7) 



^W\M^V) 



X ft 



»m(A') 



— ^ e(sLSA/S Ar )e LMAr Tr^ SL ( A0 
' u(A')=u" 

^Jm(A')' \jv{v", e) 1 ft sw(A0 k^ (A0 , V^V) 1 ) 



•L. (3.14) 



We now come to the quantization of the second integral in Eq. (13.91) . Given a graph 7 
and a 2-surface S, we can change the orientations of some edges of 7 and subdivide edges 
of 7 into two halves at an interior point if necessary, and obtain a graph 7s adapted to S 
such that the edges of 7s belong to the following four types [2]: (i) e is the up type edge 
if e fl S = 6(e) and e a (0)n a (e(0)) > where n a is the co-vector field normal to S; (ii) e is 
the down type edge if e fl S = 6(e) and e a (0)n a (e(0)) < 0; (iii) e is the inside type edge 
if e fl S = e; (iv) e is the outside type edge if e fl S = 0. In the following, we only use 
the graphs adapted to some 2-surfaces. For convenience, we will abbreviate the coordinate 
(z 1 , z 2 , !/ 3 ) of a point in 2-d surface Sy3 as (z,y 3 ). In a suitable operator-ordering, that 
integral can be quantized as an operator acting on a cylindrical function as follows: 



d 3 y 



d 2 z X W> y2 ' zl ' z2) E*{z\ z\ y 3 ) 



Ar(z,y 3 ;e') 



eeE(7) 



d 3 y 



^ Q(e',S y 3) 



[d a x 3 Ax,y)]E^y)-f, 

x 2 Ay\y 2 ; e'(0)) 



e'GE( 7 ),e'(0)GS a 3 



ir(e'(0),e') 



x [d aX s £ (x,y)] / dte a (t)6 3 (y, (e(t)))X*(t) • / 7 



1-ihKp) 2 



E lim E^( x ' e ( tfc )) ~ Xefo e (^-i))] 



n— >oo 
eeE(y) k=l 



7 



e 'es(7),e'(0)e5 e3{tfc _ i) 



/ / q x XM^fa-O.e^tfc-i); e'(O)) , , 

e(e , Aesfo-o; f ; ^ „ A e ,(0)A e (t fc _i) ■ ; 7 , 



Ar(e'(0),e') 



(3.15) 



where = to < ti < .. < t n — 1 is an arbitrary partition of the interval [0, 1], X\(t) :- 
[h e (0,t)rih e (t, l)]ABd/d[h e (0, 1)]ab (we denote X l e := X*(0) in the following) and 

-1, if e is of the up type with respect to S; 
g(e, S) = \ —1, if e is of the down type with respect to S; 

0, if e is of the inside or outside type with respect to S. 



Let us introduce the set of isolated intersection points of 7 and S 

P(y 1 S):={enS\ Q (e,S)^0,eeE( 1 )}. 



(3.16) 



The first integral in Eq. ( 13.91) can be quantized straightforwardly, since \J det(a)(x) is 
given by [4, 16] 



vSP(-y,S) V e,e'£E(j); 6(e)=6(e')=« 



veP(-y,S) 



- • 



(3.17) 



where g(e,e') := g(e, S)g(e', S). Putting Eqs. (13.141) . (13 . 1 5[) and (I3.17P together, we finally 
obtain the regularized operator corresponding to Eq. (13.91) . acting on a cylindrical function 

as 



2 9 f3 2 



81h 4 K 5 



E Ar - E E^^' 6 ^))-^^' 6 ^- 1 ))] 



wGP(7,5) eG£;(7) fe=l 
X 2^ > ^t*-i)J 7 ; A ^-e'(U)A e (t fc „i y 

Ar(e'(0), e') 



e'eB( 7 ),e'(0)£S 3 



(tfc-l) 



E ^m') 

u'ev(7) 



v(A)=v' 



x h 



h sAAy\/v(v',e) h SK{A) h s ^ (A) ,^V(v',e) ) 



x 



E 

*/'ev(7) 



K E e(^SM^)e iMJV Tr(/ isi(A , ) K^ A , y yJv{v",e) 



v{A')=v" 



8 



x ft 



sm(A') 



Now we perform the limit n — > oo, e' — > and e — * in reversed order. Keeping n fixed, 
for small enough e, only the term with k = 1 in the sum survives provided that 6(e) = v, 
and only terms with v = v' = v" contribute. So for small enough e, the above operator 
reduces to 

2 9 /5 2 



E Q,k ( S ) ' fl 



«GP(7,5) 



e'6E( 7 ), e '(0)6S e3(to) 
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^(o) £ x:(o) 

b(e)=i> 



X 



£ e(s lS js K )e IJK Tr(h 



SJ (A) 



u(A)=« 



X ft, 



aj(A) 



^WV^) 



ft.. 



:(A) 



X 



£ e(s LSA/ ^)e LMAf Tr(/ i 



^l(A') 



u(A')=u 



X ft 



ft" 1 

'V(A')' 



V>,e) 



ft 



*(A')> \M t; ' e ) 

^C(A')' V^^' 6 



A 



For small enough e', the function x^(v,e'(0)) vanishes unless w = e'(0). Hence the above 
regularized operator reduces to 



ugP(7,S) 



6(e')=u 



&(e)=i> 



e{ Sl sjs K )e IJK i:i(h Sl 



(A) 



u(A)=« 



X ft 



C(A),V^^) 



ft. 



jr(A) 



Kk{av \/V(y,e\ 



x 



<s L s M s N )e LMN TT(h SL(AI) \h^ (AI) ^V(v,e) 



v(A')=v 
X ft SM (A') 



ft 



sjv(A') 



^C(A')' \/v( V > e ) 



Notice that £>(e', 5) implies that the edges inside 5" have no contribution to the operation, 
while the action of area operator Ar(v, e') on the edges transversal to S is non- vanishing. 
Hence 1/Ar(v, e') is well defined. Thus one can take the limits and obtain an operator as 



9/32 



2 9 /3 



?. ? 

b(e)=u 



i;eP(7,5) 6(e')=« 
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Ev) 



^2 e(-S/s J s E -)e /Jif Tr^ S/(A) 



v(A)=v 



h~ l \ V 



x h 



«j(A) 



h~ x \ V 

' 4 sj(A)' V 



»jc(A) 



E(v) 



v(A')=v 
X /i SM (A') 



E <s L s M s N )e LMN Tr(h SLiAf) 



fi sz,(A')> V Vv 



sjv(A') 



'%(A')' V ^ 



A 



veP(-y,S) 



(3.19) 



Let E v ^{^) = {e G £'(7); i> = 6(e); e = * type} where * = u, d, i for 'up, down, inside' with 
respect to 5 respectively, and let X l v ^ = ^ eeB « X l- Then one can check the commutation 
relation [2] 



(3.20) 



Hence one has 



£ ^(e',5)^, KHKu-Kd)(K»+ x U+ x i,i) 

b(e')=v b(e)=v 

b(e)=v b(e')=v 



(3.21) 



Thus there is no operator-ordering problem for these two operators. Moreover, the operator 
J2 b (e')=v ^( e '' S ) X l> Eb( e )=„ x l is g au g e invariant since 



6(e\ S)Xl, X l E X e" 

b(e')=v b(e)=v b(e")=v 

[( X l,u ~ X l,d)( X l,u + X l,d + ^Qi); + X l,d + ^,i] = °- 



(3.22) 



Therefore our quasi-local energy operator Eq^(S) in (13.191) is gauge invariant. Moreover, 
we can also define a symmetric quantum version of Eq^(S) as 



(3.23) 



A special property of the Eq^S) (or Eq k {S)) is immediately clear. Because X/&(e)=i;-^e 
generates the internal gauge transformations, Eq >k (S) (or Eq k (S)) vanishes on gauge- 
invariant states. Note that the quantization of Eq ik (S) is not unique. An alternative 
quantization with a partial gauge fixing is given in Appendix [A] 
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3.2 Quasi-local normal momentum operator 

The so-called normal-directional momentum of S in adapted coordinates can be expressed 
as [20] 

Jq,i{S) ■=- [ d 2 x^det(a) I, (3.24) 
K Js 

where / is the trace of l a b, which is the extrinsic curvature tensor of S with respect to its 
unit normal u a orthogonal to E, i.e., 

I =a ab l ab = a ab V a u b . (3.25) 

Let Ar(S) be the area of a closed 2-surface S. The Euclidean Hamiltonian constraint with 
lapse function N = 1 reads 

H E (l) = ±- f d 3 x-2== e^F^El (3.26) 
Then we have (see Appendix IB1 for a proof) 

J Q , l (S) = -{H E (l),Ar(S)}. (3.27) 

Thus the quasi-local momentum Jqj(S) of S can be regarded as the "time derivative" 
of the area of S. Since there are densely defined operators corresponding to H E (1) and 
Ar(S) in Ti^, we may replace Poisson brackets by commutators times l/(ih), functions 
by operators, and obtain the operator Jqj(S). However, there are two quantum versions of 
H E (1) in the literature, the graph-changing one [19,21] and the no-graph-changing one [22]. 
It turns out that the no-graph- changing version is more convenient in defining the quantum 
version of Jq,i(S). 

We now introduce the minimal loop prescription [23]. Given a vertex v of 7 and two 
different edges e i; e,- incident at and outgoing from v, a loop within 7 starting at v 
along ej and ending at v along ej 1 is said to be minimal provided that there is no other 
loop within a%j satisfying the same restrictions with fewer edges traversed. We denote 
by L(v,ei,ej) the set of minimal loops with the data indicated. The non-graph-changing 
symmetric operator H E {1) acts on a cylindrical function as [22] 

H E C\\ f 1 ijk e ( e i> e h e fe) 

1 rJj 3iHk 2 {3 ^ ^ \L(v,e hej )\ 

r veV(-y) e i ne i ne fe=? ; 1 v ' l ' ■> M 

x ^ {{h aiV K k [h-lV v ])) ■ f, 

aijeL(v,ei,ej) 

=: &v ■ A. (3-28) 

u6V(7) 
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where {•, •} denotes the anti-commutator. 

Hence it is easy to obtain a well-defined operator corresponding to Jqj(S) as 

jQ,i(S) = ^-[H E (l)Jr(S)]. (3.29) 

It is easy to show that Jqj(S) is a gauge-invariant, diffeomorphism-covariant and symmetric 
operator. For later purposes we write a more explicit expression for Jq^S) operating on 
a cylindrical function as 

Jq,i( s ) ' a = ^ E [^?>J-/ 7 = ^ E a 

«'eV(7),«eP(7,S) «eP(7,S) 
=: E "W/y. ( 3 - 3 °) 

i)SP(7,S) 

where in the first step we used the fact that H E (1) is a non-graph-changing operator, and 
in the second step we exploited that Ar v only acts on the edges incident at v G P(j, S) 
so that the commutator with H®, , which contains only holonomies of edges incident at v', 
vanishes if v' ^ v. 



4 QLE operators 

In this section, we will use the two well-defined operators Eqj,(S) and Jqj(S) constructed 
in the last section as building blocks to quantize several types of QLE expressions. 

4.1 Brown- York energy operator 

The system under consideration is a spatial three-surface S bounded by a two-surface S 
in a spacetime region that can be decomposed as a product of a spatial three-surface and 
a real line-interval representing time. Suppose that the 2-metric a a b induced on S has 
positive scalar curvature. Then by the embedding theorem there is a unique isometric 
embedding of (S, <r a b) into the flat 3-space. Let k Q be the trace of extrinsic curvature of S 
in this embedding, which is completely determined by and is necessarily positive. The 
time evolution of the two-surface boundary S is the timelike three-surface boundary 3 B. 
Brown and York defined their QLE by the Hamiltonian-Jacobi method as [7]: 

Eby(S) f d 2 x^det(a)(k - k), (4.1) 

K Js 

where k is the trace of the extrinsic curvature k a b of S corresponding to the normal n a 
orthogonal to 3 B ( "orthogonal boundaries assumption" ) , and the integral of k Q is a reference 
term that is used to normalize the energy with respect to a reference spacetime, not 
necessarily flat. The second integral in Eq. (14.11) have been quantized as Eq. f)3.19p . The 
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construction of reference term from a reference space amounts to posing and solving an 
isomeric embedding problem. One natural choice is to embed S isometrically into Euclidean 
three space (M 3 , 5 a b) in order to obtain an extrinsic curvature tensor (k ) a b (and hence k D ). 
With the second fundamental form (k ) a b expressed in terms of the embedding's coordinate 
chart, this is the Weyl's problem, a classic problem of differential geometry for which an 
extensive literature exists [24]. However, it is very difficult to obtain the solution of (k ) a b 
in terms of a a b- Nevertheless, since the function of the reference term is to normalize the 
energy, in the quantum version one may regard it as a c-number — K = ^ J s d 2 x^/ det(a) k Q . 
In this sense the Brown- York energy has been quantized as 

E BY (S) = E Q>k (S) - K . (4.2) 

Of course, one may also take the other viewpoint that the reference term is dynamical and 
thus should be quantized. In certain symmetric models, it is indeed possible to solve k . 



4.1.1 QLE in spherically symmetric model 

We now study the Brown- York QLE in spherically symmetric quantum geometry. Consider 
a static spherically symmetric space-time with line element 

ds 2 = -N 2 dt 2 + H 2 dr 2 + R 2 (d6 2 + sin 2 6d(f) 2 ), (4.3) 

where N and H are functions of r only. Let H t be the interior of a t = constant slice with 
two-boundary S r=const specified by r = constant. A straightforward calculation gives the 
trace k of k ab as 

k = — , (4.4) 
RE ' v ; 

where the prime denotes the partial differentiation with respect to coordinate r. Now 
consider a round sphere with radius R embedded in (IR 3 , 5 a b)- Such a sphere has an extrinsic 
curvature {k ) a b with trace 

K = | ■ (4.5) 
So the Brown- York quasi-local energy can be written as 

E BY (S r=const ) = ~ d<P d6R 2 sin6 

- ^« (1 - L (4.6) 

K si 

In the spherically symmetric model, the invariant connections and triads can be written 
as [25] 

A a =A r (r)A 3 (dr) a + [Ai(r)Ai + A 2 (r)A 2 ] (d0)« + [Ai(r)A 2 - A 2 (r)A 1 ] sm6(d<p) a 
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+ A 3 cos6>(d<^) a , 
E tt =Er(r)A 3 fw9(j£j + [E\r)A 1 + E 2 (r)A 2 ]sme(^j 

+ [E 1 (r)A 2 -E 2 (r)A 1 ] (J- 



(4.7) 



where A r , A u A 2 , E r , E 1 and E 2 are real functions on an one-dimensional, radial manifold 
M with coordinate r, and Aj are the sw(2)-matrices and identical to Tj/2 or a rigid rotation 
thereof. For convenience, one introduces variables 

A v (r) := y /{A 1 {r)f + {A 2 {r))\ (4.8) 
E*{r) := ^{E\r)f + {E\r)f, (4.9) 

and a(r), f3(r) defined by 

A lC os/3(r) +A 2 sin/5(r) = (^i(r)A 2 - 4 3 (r)Ai) /A v (r), (4.10) 
Ax cos (a(r) + f3(r)) + A 2 sin (a(r) + f3(r)) = (E 1 (r)A 2 - E 2 (r)A 1 ) /E v {r). (4.11) 

In the spherical coordinate system, the components of the spatial 3-metric qw, on E take 
the form 



diag 



(^£,^,^sin^). 



So we have the relation 



The QLE reads 



H 



E v 
\/W T 



R 



'E r . 



E B y{S; 



r=const ) 



8vr 



E T 



2E^ 



(4.12) 



(4.13) 



(4.14) 



In Ref. [26], a canonical transformation from (A r , A u A 2 ; E r , E 1 , E 2 ) to (A r , /3K V , rj; E r , E* ', P*) 
where 



K v {r) := v / (^i(r)) 2 + (K 2 (r)) 2 

is the extrinsic curvature component and 

rf(r) := a(r) + f3(r), 
i">(r) := 2A tp (r)E' fi (r) sin a(r) 

has been made. The symplectic structure is then given by 

Sl B = % I dr{SA r A SE r + 2(35K V A bE v + 5r] A 5P r ') 
K P Jb 
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(4.15) 



(4.16) 
(4.17) 



(4.18) 



The resulted Hilbert space is spanned by an orthonormal basis of spin network states: 

T g ,k,,x(A) = JJexp ( \ik e jA r (r)dr J JJ exp(ifi v f3K v (v)) exp(ik v r](v)) (4.19) 

with edge labels k e G Z, vertex labels G K. and fe„6Z for graphs g in the 1-dimensional 
radial manifold M. The flux operators corresponding to the momenta E r and E v act on 
the spin network states as: 

gH^^ ^a r rt> , (4.20) 



f Hr)T g , K „ = ^ ^^(r, !,)^, (4.21) 



where e ± (r) are the two edges (or two parts of a single edges) meeting at r, 8(r, v) is the 
Dirac delta distribution, and X is a region of the reduced (radial) manifold M. 

We now consider the quantization of the QLE. The spin connection component can be 
regularized as [26] 



2^ 4 I P+ Mr f E^dr 

Now the E r in Eq. (021) can be promoted as well-defined operators, and the j E^dr ma y 
also become well-defined operator by suitable treatments [27,28]. Hence the quasi-local 
energy operator Eby (S r = CO iist) in the spherically symmetric sector of LQG can be well 
defined. 

4.2 Liu-Yau energy operator 

Let I be the trace of the extrinsic curvatures of the 2-surface S in the physical spacetime 
corresponding to the future pointing timelike normal. Liu and Yau define a quasi-local 
energy by [8] 

E LY {S) :=- [ d 2 x.ffet{^)(k - J\k 2 -P\). (4.23) 
« Js 

Since the first term in Eq. (14.231) is again a reference term, we will only consider the second 
term 

Eqm{S) = -- [ d 2 xy/^t{a)^\k 2 -P\. (4-24) 
K Js 

Let g e (x, y) be a 1-parameter family of fields on S which tend to S 2 (x, y) as e tends to zero, 
i.e., such that 



lim / d 2 yg e (x 1 ,x 2 ;y\y 2 )f(y 1 ,y 2 ) = f(x\x 2 ) (4.25) 
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for all smooth densities / of weight 1 and of compact support on S. (Thus, g e (x,y) is a 
density of weight 1 in x and a function in y.) Using g e (x, y) as smearing function, one can 
regularize E Q:k j(S) as 



d x 



d 2 x 



- / d 2 yg e (x,y)y/det(a)(-k + I) - / d 2 zg e (x, z)y/det(a)(-k - I) 



1/2 



[E Q - k+l ] g ^x) [£ Q) _ fc _,]^( 



1/2 



d 2 xJ[EsUx). (4.26) 



It is easy to see that En kl (S) tends to EQ jkj i(S) as e tends to zero. Let us now turn to the 
integrand of Eq. (I4.26p . It can be promoted as an operator acting on a cylindrical function 

as 



[ E s]g e { X ) ■ fl = \[ E Q-k+l\gX x )[ E Q-k~l]gX X ) 



= 9e(x 1 v)g,(x,v')\(E QAv + Jq,i, v )(E QAv > - J Q ,iy)\ ■ f 7 , (4.27) 

where the absolute value \(EQ ik , v + Jq,i,v){Eq^ )V — Jq,i,v)\ = A indicates that one is sup- 
posed to take the square root of the operator A^A. We choose e sufficiently small so that 
g e (x,v)g e (x,v') is zero unless v = v' . Then one obtains 



i E s]gS X )-fj= X] 9e(x,v) 2 \(E QAv + J Q ,l jV )(E Q>ktV - J Q ,l lV )\ ■ f T 



(4.28) 



Notice that [Es} g ,,(x) is a non-negative self-adjoint operator and hence have a well defined 
square root, which is also a positive-definite self-adjoint operator. Since we have chosen 
e to be sufficiently small, for any given point x in S, g e (x,v) is non-zero for at most one 
vertex v. We can therefore take the sum over v outside the square root and thus obtain 



[ E s]gS X ) ■ A = 9e(x,v)\/\(EQ Av + J Q ,l :V )(E QAv - J Q ,l, v ) \ • / 7 - 

veP(-y,S) 



(4.29) 



Finally, we can remove the regulator. By integrating both sides of Eq. (I4.29P on S and 
then taking the limit e — > 0, we obtain the desired operator corresponding to Eq. (I4.24p as 



E Q,k,l(S) ' fj — ~ \/ (EQ,k,v + JQ,l,v)( E Q,k,v - Jq.I,v) 



(4.30) 



4.3 Hawking energy operator 

By studying the perturbation of the dust-filled open Friedmann-Robertson- Walker space- 
time, Hawking found that [9] 



1 



1 



167T 



d 2 xy/det(a) (k 2 - I 2 ) 



(4.31) 
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behaves as an appropriate notion of energy surrounded by the space-like topological 2- 
sphere S. The virtue of Hawking energy is that it does not need a reference term. We can 
regularize it as 



E H (S) =^^/MSj- lim y^^l f d 2 x^/d^)(-k + 

K e-*Q 8y/lTK Jg 



I) 



x / d 2 y- 



s e 



jdMj(-k-i) 

! V / detR V 



: V^^(S)-lim^5^ / *xy/tet(a)(-k + I) 

K e-»0 8V7TK J s 



x / d 2 y 



Ar(x, e 



(4.32) 



To quantize the expression (14.321) . one can replace the yj Ar(S) by y Ar(S) and use the well- 
defined operators Eq^ and Jq,i, v - We then formally get an operator acts on cylindrical 
functions as 



,l,v 



Ar{S) 
Ar, 



■x 



veP(rf,s) 

x X] u ')(^Q,fcy - Jq,i,v>) ■ f- 



7' 



(4.33) 



i/eP(7,s) 

For sufficiently small e, ( w ? v is zero unless v = v' . Thus one has 

E H {S) ■ /y 



2^F 



Ar(5 • / 7 - — =r- V £ gfc +j g , V (E Qtk , v -J Q>ltV )-f 7 . (4.34) 



ueP(7,s) 



However, this is not a densely defined operator due to Fortunately, since Eq 



k.v 



vanishes the internal gauge-invariant states, for a gauge-invariant state \l/ 7 e H D we have 



2^ 



K 



(4.35) 



^GP( 7 ,S) 

Hence we can obtain a well-defined operator corresponding to Hawking energy in 7i as 
E H (S)-%=^JAr(S)-^ 



+ 



H^yJAr v 



Ar(S) 



It is easy to see that this operator is symmetric in 7i Q . 



H:,\^Ar v 



# r (4.36) 
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4.4 Geroch energy operator 

Geroch modified the Hawking energy and gave the other definition for QLE as [10] 

1 



E G (S) := 
It can be regularized as 



Ar(S) 
16ttG 2 



1 



16vr 



d 2 x^/det(a)^ 



(4.37) 



Eg(S) =-^—JAr(S) - km v ; / rf 2 z v , . v / k / d' 
The quantum operator corresponding to the Geroch energy then formally reads 



yx 2 e (x,y)y/det(a)k. (4.38) 



E G (S) ■ U A) 



K 



Ar{S) ■ UA) 



^T\/^) E j-^-fM)- (4-39) 



^GP(7,5) 

Let T s be the gauge-invariant spin network function in H . Then Eq(S) acts on T s as 

2 V? 



H 



K 



Ar v -T s . 

veP(j,S) 



(4.40) 



Thus Eg{S) is well defined in TL Q . Hence the spectrum of the area operator [16] implies 
the spectrum of Eq{S) as 



Spec[£ G (S)] 



2nh(3 



E + 1) + 2jr or + 1) - jr uj o^ 

veP(%S) 



(d+u)/ -(d+ii) 



where and ji d+u ^ are half- integers subject to the usual condition: 



■(d+u) _ r | -(d) _ ■(«)! I -(d) _ ■(«)! , i „•(<*) , Au)\ 

JV C \ W8 Jl> I ) I Jl! Jl> T ) ' ' ' )Jt) ^ Jll J • 



;(<*) 



;(«)! U( d ) 



Ml 



;(<*) 



(4.41) 



(4.42) 



In the case that all edges of the graph 7 underlying T s puncture S, i.e., 7 has no edges 
tangential to S, the spectrum is reduced to 



Spec[E G (S)] 



E y/jvtiv + 1) 

ueP(7,s) 



E a/jUj'i- + K 

veP(y,S) 



(4.43) 



where j v are half-integers. In summary, Eg{S) is a densely defined, positive semi-definite 
operator in Ti. Q , and its spectrum is entirely discrete inherited from the property of area 
operator. Thus we have proved a quantum positivity QLE theorem. Moreover, Eq(S) 
is both internal gauge invariant and invariant under the diffeomorphism transformations 
tangent to S. Furthermore, the Geroch gravitational energy in a quantum state labelled 
by a graph 7 is concentrated at the vertices of 7 which live on 2-surface S and the edges 
which puncture S transversely. The discreteness of quantum gravitational energy enable 
us to estimate the statistical entropy of the region enclosed by S in next section. 
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5 Discussions: entropy-area relation in LQG 



It was first speculated by Bekenstein that one could associate an entropy &bh to a black 
hole with horizon area A as [29] 

where ks is the Boltzmann constant. The statistically mechanical origin of this entropy 
has been an outstanding mystery for physicists. Some intuitive arguments and accurate 
calculations have been done in the framework of LQG to account for Eq. (15. If) [30-33]. 
Now we have the quantum gravitational energy of any finite region bounded by a closed 2- 
surface S. So, in principle, one may study the thermodynamical properties of an arbitrary 
bounded gravitational system in LQG by the standard statistical mechanics method. For 
example, if one considers it as a canonical system, the partition function reads 

Z(S) = Tr e"^ (5.2) 

for certain QLE operator E(S). Then all the thermodynamical quantities including the 
entropy of the system can be derived in principle. However, since the spectrum of the 
above QLE operators are either too complicated or unknown yet, it is still difficult to do 
practical calculations. Thus, further investigations in this strict approach are needed to 
understand the spectrum properties of the QLE operators. 

On the other hand, the so-called holographic principle says that, at the fundamental 
(quantum) level, one should be able to characterize the state of any physical system lo- 
cated in a bounded spatial domain by degrees of freedom on the surface of the domain. 
Consequently, the number of physical degrees of freedom in the domain is bounded from 
above by the area of the boundary of the domain instead of its volume. If the entropy 
representing the degrees of freedom including gravity in a bounded domain could be cal- 
culated, one would be able to check whether the holographic principle is valid or not in 
LQG. Some QLE operator and the corresponding partition function provide a possible 
approach to this issue. Again, we need more control on the spectrum of the QLE operators 
constructed above. 

Nevertheless, the virtue of Geroch QLE operator Eq(S) in (14.401) is beneficial for us to 
generalize the entropy-area relation in the framework of LQG. Our discussion is restricted 
a simple self-gravitating system bounded by a space-like closed 2-dimensional surface 5". 
One can count the number of quantum states corresponding to eigenvalue of E G (S). Thus, 
it is regarded as a microcanonical ensemble where the energy of the system is fixed. The 
spectrum ( I4.43[) of the Geroch energy operator implies that the energy eigenvalue involves 
only the number N of punctures and the spins j v of the edges that intersect the surface 
S. Thus the number of the eignstates of a given eignvalue of Eg(S) is infinite because 
different positions of the punctures give different states. This is no longer the case after 
modeling out the spatial diffeomeorphisms tangent to 5*. Following [30,34], we shall treat 
the punctures as distinguishable. Our task is to count the number M(Eg(S)) of quantum 
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states corresponding to the classical QLE Eq(S). It is to see that Af(Ec(S)) is the same 
as the number Af(Ar(S)) corresponding to the classical area Ar(S) = f^E G (S) 2 . Hence 
the states which we are considering satisfy 

Ar(S) = Kh/3^n jy /j{jTlj, (5.3) 

3 

where rij is the number of punctures with spin j. Following [33], the number of the states 
is given by 

Af(Ar(S)) = n(2j + IP ■ (5-4) 

[[ j" r 3 

Using Stirling's formula, one gets the entropy as [33] 

§ = k B \nAf(E G (S)) = k B ]nJ\f(Ar(S)) = k B ^^j , (5.5) 

where /3 is the solution of the equation 

1 = ^(2j + l) e - 2 *P°V3W~) . (5.6) 

3 

Moreover, if one assumes that N(Eg{S)) represent the physical degrees of freedom in 
the domain bounded by S, the holographic principle is realized in LQG. Note that for 
marginally trapped surfaces, the QLE must be the irreducible mass tJA(S)/16ttG 2 [1]. 
Hence our discussion for these cases do not depend on a specific definition of QLE. 

We conclude with a few open issues related to the present work, (i) More understanding 
of the spectrum of the QLE operators is needed in order to do further practical calculations, 
(ii) Semiclassical analysis on the QLE operators is yet to be done, (iii) Since classically the 
integral of k is non-zero in general, the vanishing of QLE-like operator Eq^(S) on gauge 
invariant states is in some sense awkward. One may consider other possible operator 
orderings in its construction in order to avoid the weakness, (iv) The knowledge on the 
physical Hilbert space of LQG would be a great help to our scheme. 
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Appendix 



A The partial gauge fixing version of Eq^(S) 

Classically, we have 

k ^ = D ° {-7m n ' K ) = 7m aMEt) ' (A1) 

where we have introduced the internal vector n l = n a ef in the "internal space". It is 
convenient to first carry out a partial gauge fixing. Let us fix an internal vector field n l 
with n l rii = 1. We restrict ourselves to flat derivative operators d which annihilate n l in 
addition to 5ij. We call the partial gauge fixing as n l -gauge fixing. Fixed n gives a fixed 
direction in the internal gauge group SU(2). Thus the structure group is reduced from 
577(2) to U(l). Physical states and observables should, of course, be independent of this 
choice. The n l gauge transformations on the 2-d surface S, which keep n % invariant, are 
generated by the following Gauss constraint 

G{\n l ) = J d^xXn'Gi, (A.2) 

where A is an arbitrary real number. Hence the corresponding Gauss constraint operator 
is given by 

Gf(An') • / 7 = -^An' £ X\ • / 7 . (A.3) 

ee£(7) 

Let ifi^ be a n*-gauge-invariant cylindrical function corresponding the above Gauss con- 
straint on A. Then at every vertex v of 7, the following condition must hold: 

n* X i-^ = °- ( A -4) 

where X l v = J2 e eE(-y) We)=v Under the n*-gauge fixing, the extrinsic scalar curvature k 
reduces to 

d a Ef. (A.5) 



Note that now n l is a non-dynamical constant which need not to be quantized. For sim- 
plicity, we choose again adapted coordinates {x 1 ,^ 2 ,^ 3 } with respect to S. Our aim is to 
quantize the following quantity under the partial gauge fixing 



E Qk (S) = - 1 - [ ^ V^gg(^) n\ X )d a E?{y) 
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— lim - 

e^O K J S 



det(cr)(x) 



lim — 



d X 



det(q)(x) 
y/det(a)(x) 



d 2 x v ^. K > n\x) I A^i^ d a E?(y) 



V(x,e) 



E 



n\x) / d 6 y{d aX t{x,y)}Et{y) 



(A.6) 



The first integration of Eq. (1A.6I) can be written as 



d xn l (x) 



y/det(a)(x) 
V(x,e) 



d xn l (x)> 



'n a (x)Ef(x) n b (x)Ef(x) 



V(x,e) 



(A.7) 



where n a (x) = (dx 3 ) a . Using g e '(x,y) satisfying Eq. (14.251) as smearing function, we define 



~n a (x)E?(x)~ 




[ V(x,e) 


9 t i ^ S 



2 , n a (u)EHu) 
d ug e i{x, u) — - 



V(u,e) 



S V l M 7 € ) L 



f<- 2 I 



e ,(x,u)e ijk {Ai(u), ^V(u,e)} A {A h (u), ^/V(u,e)}, (A.8) 



where we have used the classical identity 

ei(u) = -{A i a (u),V{u,e)}, 

K 

and absorbed 1/V(u,e) into the Poisson brackets. Thus we have 
y/det(a)(x) 



d xn l (x)- 



V(x,e) 



e'-> 



lim / d 2 xn l (x 



~n a (x)Ef(x) 




h b {x)E\{x) 




[ V(x,e) 


9 e > 


[ V(x,e) 


9 f > 



(A.9) 



(A.10) 



We introduce a triangulation of the 2-d surface S [35,36]. Denote by A a solid triangle. 
Single out one of the corners of the triangle and call it v(A). At v(A) there are incident 
two edges S\(A), 52(A) of dA which we equip with outgoing orientation, that is, they start 
at v (A). Let us now write the integral over S as a sum of integrals over A where A are 
triangles of some triangulation T of S, 



n a (x)E?(x) 
V(x,e) 



J 9 e i 



-J>(*XA)) £ 



e(s I sj)e IJ e ijk x 



Sj (A)n Sj (A)=„(A) 

Tr (^^(A)^- 1 ^, y/V(v(A),e) }) Tr (r k h Sj(A) {h^ {A): y/V(v(A),e) } 

(AH) 

where e(sjSj) := sgn(det(i/Sj)(w(A))) takes the values +1, —1, if the tangents of the two 
segments sj, sj at v (in that sequence) form a matrix of positive, negative or vanishing 
determinant. For convenience, we introduce 



(A.12) 



4(v(A)) := Tr (r 3 h Sl{A) {h^ {AV y/V(v(A),e) }) . 
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Then the quantity in the square root of the expression in (lA.lOj) can be regularized as 



~h a (x)E?(xy 




n b (x)E h i (x) 


[ V(x,e) 


9e' 


[ V(x,e) \ 



9e' 



1 6 

= rf E ^(^^(A))^(^^(A , ))e(^^)^e / V / (t;(A))e^(^(A)) 

A,A'eT 

x e(s K s L )e mn e KL e%(v(A'))e n L (v(A')). (A.13) 

We introduce a triangulation of the 2-d surface 5* in adaption to the graph 7 [35, 36] 
and only consider the terms in flA.13j) which sum over triangles A whose basepoint v (A) 
coincides with a vertex v of the graph, 

^ 5^ 9e'(x,v)g e ,(x,v')e ijk e imn ^ -^—e(s I sj)e IJ e 3 I (v)e k J (v) 

x E ^e^^e^'KK)- (A.14) 

sx(A')n Si (A')=^ ^ ' 

For sufficiently small e', g e i(x,v)g e '(x,v') is zero unless t> = v'. Then the double sum over 
vertices reduces to a single one, Eq. (1A.14D reduces to 

16 ^— \ r 1 m2 16 

— 2^ i9,'{x,v)\ e ljk €i 



«ey(7)n5 v ' 

sj(A)nsj(A)=Sic(A)nsi(A)=u 

= f E [&'M] 2 ;J^ E e(^^)e(^^)e /J e X M( W )e^( W )e^( W )ei(tO. (A.15) 

where we have abbreviated E S/ (A)n Sj (A)= SK (A)n Si (A)=t, -> J2ijkl- We further introduce 
the manifestly gauge invariant quantities [37] 

Qik{v) = ^{vy K {v). (A.16) 

The (1A.14I) can be written as 

I £ [^(x,i;)] 2 -^- E ^/^M^^W^jl^)- (A.17) 

To write the regulated operator corresponding to (1A. 13[) . we replace Poisson brackets by 
commutators times l/(ih), holonomies by multiplication and V by V. As we evaluate the 
operator corresponding to Eq. flA.13j) . we find out that only those triangles A contribute 
whose basepoint v(A) coincides with a vertex v of the graph due to the presence of the 



23 



volume operators. Hence we obtain the regulated operator corresponding Eq. (1A.13j) acting 
on a cylindrical function as 



~n a (x)E?(x)~ 




n b (x)E\(x) 


[ V(x,e) 


9 t > 


[ V(x,e) \ 



f- 



32 1 



16 

Yl l9e'(x,v)] 2 —— Y t(siSj)e(s K s L )e IJ e KL q IK (v)q JL (v) ■ f 1 



(my 

^2 [g e >(x,v)] 2 q v • /. 



E(v 



IJKL 



(A.18) 



where 



with 



q IK (v) = e\{v)e 3 K {v) 



= Tr I rj/i Sj(A ) 



uGV(7) 



(A.19) 



(A.20) 



It is easy to see that the operator in flA.18j) is gauge invariant. Notice that the self- 
adjointness of ie](v) implies that q v is a non-negative self-adjoint operator and hence has 
a well defined square root. Since we have chosen e' to be sufficiently small, for any given 
point x in S, g e '(x, v) is non-zero for at most one vertex v. We can therefore take the sum 
over v outside the square root and obtain 



n a {x)Ef{x) 
V(x,e) 



9e' 



h b {x)E\{x) 
V(x,e) 



1/2 



9t' 



A = 9e>{x,v)yfq~ v ■ / 7 . (A.21) 

vev(-y)ns 



Finally, we can remove the regulator, i.e., take the limit as e' tends to zero. Then the 
following equality holds in the distributional sense. 



y/det(a)(x) 
V(x,e) 



v&V{~ f )C\S 



The second integration of Eq. (1A.6I) can be similarly quantized as 



d 3 y[d aX 3 e (x,y)}E?(y) • / 7 



2 ' ' n— >oo 

eG£(7) k=l 



(A.22) 



V lim V M(x,e(4))-Xe(^e(Vi))l^M-/ r (A.23) 
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Using Eqs. (1A.22I) and (1A.23I) . we obtain the regularized operator corresponding to Eq^(S) 
as 

^ s vev(-y)ns 

n 

x E E[^( x ' -^( x ' e ^-0)] x e(vi)-/ 7 

eeE(7) fe=l 



2 

v£V(y)nS 



n (v) 



n 

X 

eG E( 7 ) fe=l 



n 

E E [Xe(v,e(t k )) - xf^.e^-i))]^^!) • / r (A.24) 



Now we perform the limit n — > oo and e — > in reversed order. Keeping n fixed, for small 
enough e only the term with k — 1 in the sum survives provided that 6(e) = v. We then 
obtain the operator 

e q AS) • A = ^ E [^M") E x e • A 

t)GV(7)n5 6(e)=u 

= £ &] 1/2 n>) E y e-/7, (A.25) 

«GV(7)n5 b(e)=v 

where Y* = —\X % e is the self-adjoint right-invariant vector field. It is clear that the Eq^(S) 
in (IA.25I) is densely defined in 7ikm, and it vanishes the n*-gauge-invariant states. 

B Proof of an equality 

We first give a proof for Eq. (I3.27P in 13.21 By choosing adapted coordinates {x 1 ,^ 2 ,^ 3 } 
with respect to S, the normal-directional momentum of S is given by 

JqAS) :=- / d 2 xy/d^)l 
K J s 

=- / d 2 xJh c n c det(g) (K - n a n b K ab ) 
k Js 

=J 1 (S) + J 2 (S), (B.l) 

where 

/ :=a ab l ab = a cd V c u d = V a u a - n c n d V c u d . 
The first term in Eq. (1B.1I) can be written as 

J^S) = - I d 2 x^/^E°;Kl (B.2) 
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and the second term as 



JJS) = - - I d 2 xJf^ c n a n b K a E b . 
« Js 



Thus we obtain 



Jq,i(S) =-( f (Pxy/fitfPKiE? - [ d 2 xy^n a n b K a E b 
k^Js Js 



On the other hand, 



Note that the area of S can be written as 

Ar(S)= J d 2 x^n a n b E a i E b b i i = J d 3 x^Jh a h b E?E b 5 i i S(x 3 , 0). 



One gets 



SArjS) n b E^h a 5 a d 5l ti b e bl n d ^ 3 



5H E {l) l 



(e abd d a e bl + e abd e i:jl Ai4) 



5A l d (y) k 

Plugging Eqs. flB.7)) and (1B.8I) into (IB.5I) . we obtain 

{H E (l),Ar(S)}= f d'y-^ie^daeu + e^A^n^hJiy^O). 

The two terms in (lB.9j) can be reduced respectively to 

e abd d a e bl h c e% =e^ d (-T^ ab e fl - e lmn T™e n h )h c e l c fi d = -T™h c h d e abd e fbc 
= - T™E a m h d fi d + T™E d m n d n\ 
e a6c e^A^n c e'n d =A\n c fi d e abc e fbc e{ = A™E a m h d n d - A™E d m fi d n a . 
Hence one has 

{H E (l),Ar(S)}= [ d*y-^{K™E a m n d h d -KEim a ¥{y\U) 



d 2 yKTE« m V^~ d - / A Q 

Js Vn e n e 

The second integrand in fIB.llj) reduces to 

K « E ^f =^n«n b KlE b . 
\/n e n e 

We thus obtain 

{H E (l),Ar(S)}= [ d 2 y^fK b K a Et- [ d 2 Vv ^ n a n b K a E b 
Js Js 

Comparing (1B.4j) to (IB. 121) we complete the proof. 
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